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At finite temperature the distribution of the total momentum is an observable characterizing the 
thermal state of a field theory, and its cumulants are related to thermodynamic potentials. In 
a relativistic system at zero chemical potential, for instance, the thermal variance of the total 
momentum is a direct measure of the entropy. We relate the generating function of the cumulants 
to the ratio of a path integral with properly shifted boundary conditions in the compact direction 
over the ordinary partition function. In this form it is well suited for Monte-Carlo evaluation, and 
the cumulants can be extracted straightforwardly. We test the method in the SU(3) Yang-Mills 
theory, and obtain the entropy density at three different temperatures. 
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1. Introduction 

Thermal field theory is the theoretical tool for computing properties of matter at high temper- 
atures and densities from first principles. It allows one, for instance, to determine the equation of 
state of Quantum Chromodynamics, which in turn is an essential ingredient to understand the prop- 
erties of matter created in heavy ion collisions, and to model the behavior of hot matter in the early 
universe (for a review at this conference see Ref. [p]]). Obtaining first-principles predictions from 
a thermal field theory is often challenging since it describes an infinite number of degrees of free- 
dom subject to both quantum and thermal fluctuations. Even though there are several established 
methods to compute the thermal properties of field theories ^ ||, |], Q], new theoretical concepts 
and more efficient computational techniques are still needed in many contexts particularly when 
weak-coupling methods are inapplicable. 

Recently Meyer and I proposed a new way to determine the equation of state of a thermal 
field theory, and it is the aim of this talk to discuss the results obtained in these papers []6[ [7|]. We 
related the generating function of the cumulants of the total momentum distribution to a path inte- 
gral with properly chosen shifted boundary conditions in the compact direction normalized to the 
ordinary one. By exploiting the Ward Identities (Wis) associated to the space-time invariances of 
the continuum theory, the cumulants can be related in a simple manner to thermodynamic poten- 
tials. In a relativistic theory at zero chemical potential, for instance, the variance of the momentum 
probability distribution measures the entropy of the system. 

Crucially the argument holds, up to harmless finite-size and discretization effects, also in a 
lattice box. The formulas are thus applicable at finite lattice spacing and volumes, where ratios of 
path integrals can be determined by ab initio Monte Carlo computations. As a result the entropy 
density, the pressure and the specific heat can be obtained by studying the response of the system 
to the shift, and no additive or multiplicative ultraviolet-divergent renormalizations are needed in 
taking the continuum limit. We have numerically tested our proposal in the SU(3) Yang-Mills 
theory, for which the entropy density has been determined at three different temperatures. 

2. Momentum distribution from shifted boundaries 

For an Euclidean field theory at a temperature T = 1 /Lq, where Lo is the length of the "time" 
direction, the relative contribution to the partition function of states with total momentum p is given 
by 

gOg) = Tr{ g - L ° H P(P)} 
V Tr{e- L <>H} 

where V = L? is the spatial volume, and L is the linear dimension in the three spatial directions. 
The trace in Eq. ( |2.1[ ) is over all the states of the Hilbert space, p( p ) is the projector onto those 
states with total momentum p, and H is the Hamiltonian of the theory. If we introduce the partition 
function 

Z(z) = Tr{e- iof V'P z } (2.2) 

in which states of momentum p are weighted by a phase e' p z , and we use the standard group theory 
machinery (see Refs. ||8|, |9|] for a detailed discussion of this point) it easy to show that 

*(p) = \\ dhe-**Z(z) . (2.3) 
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The generating function K(z) of the cumulants is defined as 



^ (Z) = ^I^(P)> (2-4) 



and the cumulants are given by 

d Zm d 2n 2 d 2n, 



k 



v 



{2n,,2n 2 ,2« 3 } V ' ^2n\ ^2n 2 g^2n 3 y 



(2.5) 

z=0 



where they have been normalized so as to have a finite limit when V — > °°. The shifted partition 
function Z(z) can be expressed as an Euclidean path integral with the field satisfying the boundary 
conditions 

0(Lo,x) = ±0(O,x + z) (2.6) 



with the + (— ) sign for bosonic (fermionic) fields respectively. From Eqs. ( |2.3D and ( |2.4| ), the 
generating function can thus be written as the ratio of partition functions 

= ^ , (2.7) 

i.e. two path integrals with the same action but different boundary conditions, and the cumu- 
lants can be obtained by deriving it with respect to the shift parameter z an appropriate number of 
times. Being the cumulants the connected correlation functions of the charges associated to the 
translational invariance of the theory, the functional K(z) and the momentum distribution R(p) are 
ultraviolet finite, see Refs. []6[ [7]] for a comprehensive discussion of this point. 

2.1 Extension to the lattice 

When defined on a lattice, the theory is invariant under a discrete subgroup of translations and 
rotations only, the momenta are quantized, and the continuum Wis are broken by discretization 
effects. Generic lattice definitions of the energy-momentum tensor, as well as the corresponding 
charges, require ultraviolet renormalization. It is still possible, however, to factorize the Hilbert 
space of the lattice theory in sectors with definite conserved total momentum. The formula for the 
lattice projector is given by 

n 3 

p(P) = —Y^ e •^ P 'V' pz (2.8) 



z 



where the sum is over all the lattice points. Since only physical states contribute to the symme- 
try constrained path integrals in Eq. ( |Q| ) when Lq ^ 0, the lattice momentum distribution R(p) is 
expected to converge to the continuum universal one without the need for any ultra-violet renor- 



malization. The definition of the cumulants in Eq. (2.5) is thus applicable at finite lattice spacing, 
provided the derivatives are replaced with their discrete counterpart, and no additive or multiplica- 
tive ultraviolet-divergent renormalization is needed for taking the continuum limit. 

2.2 Extension to other symmetries 

The factorization of the Hilbert space into sectors whose states have definite transformation 
properties under a symmetry of the lattice theory is more generally applicable than just for trans- 
lations. For a generic discrete group, the cubic rotations for instance, the projector onto the states 
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which transform as an irreducible representation /I is given by 



p(M) 



'fix 

o 1=1 



i)f(«i). f(^-)i0) = ir'}, 



where is the dimension of the representation, g the order of the group, X i^i) tne character of 
the irreducible representation for the group element and T(Mi) is the representation of the group 
element onto the Hilbert space. The generalization to a continuum symmetry, such as for instance 
the baryonic number, is straightforward. The connected correlation functions of the corresponding 
charges can be extracted from the symmetry constrained path integrals, defined analogously to 
Eq. (2.1), by studying the response of the system to the properly chosen twist in the boundary 
conditions 

3. Continuum Ward identities and connection to thermodynamics 

In the continuum theory the invariance under space-time translations implies the Wis 



S v (d^ v (x) Ox . . .O n ) = - £ (01 . . . S*Oi . . . O n ) 

t=l 



(3.1) 



where O, is a generic local field, 8*Oj is its variation under the local transformation parameterized 
by £ v (z) = E v 8^'(z — x), and T^ v is the energy-momentum tensor (see Ref. [|7|] for unexplained 
notation). By choosing e v = 8 vk e k (no summation over k), Tqq{x 1 ) as interpolating operator, the 
WI in Eq. ( 3T ) and translational invariance lead to 

d£ {f 0Q (x l )f 0k (x 2 )) c = df (f kk (x 2 )7w(x l )) c , f^ v (x) = ( \ n dx p 



T^ v (x) . (3.2) 



By choosing e v = 5 v o£o, Tok(x 2 ) as interpolating operator, and thanks to translational invariance 
and the symmetry of To k , the WI in Eq. (3T ) gives 

d£ (%o(x l ) To*(x 2 )) c = djf (f ok (x l ) f Qk (x 2 )) c . (3.3) 
By putting Eqs. (|0|) and ( ^^ ) together we arrive to 

d£{(f 0k {x 2 )f Qk (x l )) c - {f kk (x 2 )f m {x l )) c } = . (3.4) 
By integrating in x\, while keeping all insertions at a physical distance (x' all different), we obtain 

{(To fc (xJ)rofe(4))c-(roo(xi)T te (4)} c } =l\ {(7 0k (x l )7 Qk (x 2 )) c -(7 0Q (x l )f kk (x 2 )) c } , (3.5) 

where To k (xo) = j d 3 xTo k (x). If we remember that in the Euclidean the momentum operator maps 
to p k — > —iTo k , the pressure maps to p = (T kk ), and if we note that the r.h.s. of Eq. (3^) vanishes 
in the thermodynamic limit, we arrive to infinite-volume relation 



k{o,o,2} =T2 gfP 



(3.6) 
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Combined with the infinite- volume WI s = -gfP, where s is the entropy density, it leads to [||, |7|] 



1 



1 d 



By following an analogous derivation, it is possible to show that the specific heat is given by 

Id 4 3d 2 



(3.7) 



c v = lim - 



3T 4 dz 4 T 2 dz 2 



lnZ({0,0,z}) 



z=0 



(3.8) 



The last two equations make clear that the response of the partition function to the shift z is gov- 
erned by basic thermodynamic properties of the system, and that the potentials entering the equa- 
tion of state of the thermal theory can be extracted by rather simple formulas. Crucially the con- 
vergence to the thermodynamic limit is exponential in ML, where M is the lightest screening mass 
of the theory 0. 

4. Numerical computation 

We have tested the numerical feasibility of the computational strategy presented above in the 
SU(3) Yang-Mills theory. The lattice theory is set up on a finite four-dimensional lattice with a 
spacing a and periodic boundary conditions in the space directions. It is discretized by the standard 
plaquette Wilson action 



S[U] 



II 

x \i<v 



1 

1 - -ReTr 

3 



{^vW} 



(4.1) 



where the trace is over the color index, and go is the bare coupling constant. The plaquette U ilv {x) 
and the path integral Z are defined as usual. The most straightforward way for computing the 
cumulant generator is to rewrite it as 



Z fJiF(z,r m ) ' 



(4.2) 



where a set of (n + 1) systems is designed so that the relevant phase spaces of successive path 
integrals overlap, and that iF(z, ro) = Z(z) and 3f(z, r n )=Z. The path integrals of the interpolating 
systems that we have implemented are defined as 



Z(z,r) = DUDU MLo/ayi e 



-S[U,U 4 ,r] 



(4.3) 



where U^^ Lo / a ^_i is an extra (5 th ) temporal link assigned to each point of the ((Lo/a) — 1) time- 
slice. The action of each system is 



S[U, U 4 ,r]=S[U] + tL(l-r)Y, ReTr { U 0k ((Lo/a) 

3 x,k 1 



l,x)-U 4k ((L /a)-l,x)} , (4.4) 



with the extra space-time plaquette given by 
U 4k ((Lo/a)-l,x)=U 4 ((Lo/a)-l,x)U k (0,x + z)ul((Lo/a)-\,x+k)ul((Lo/a)-\,x) . (4.5) 
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Lat 6/gl L /a L/a \{^f r /a K(z) 



A! 


5.9 


4 


12 


2 


4.48(5) 


17.20(11) 


5.10(3) 


Al a 


5.9 


4 


16 


2 


4.48(5) 


40.71(15) 


5.089(19) 


A 2 


6.024 


5 


16 


2 


5.58(6) 


13.05(10) 


4.98(4) 


A 3 


6.137 


6 


18 


3 


6.69(7) 


7.32(8) 


4.88(6) 


A 4 


6.337 


8 


24 


4 


8.96(9) 


4.32(16) 


5.12(19) 


A 5 


6.507 


10 


30 


5 


11.29(11) 


2.62(17) 


4.9(3) 


Bi 


6.572 


4 


12 


2 


12.28(12) 


22.22(11) 


6.58(3) 




6.572 


4 


16 


2 


12.28(12) 


53.47(16) 


6.684(20) 


B 2 


6.747 


5 


16 


2 


15.34(15) 


17.11(15) 


6.53(6) 


B 3 


6.883 


6 


18 


3 


18.14(18) 


9.61(9) 


6.40(6) 


B 4 


7.135 


8 


24 


4 


24.5(3) 


5.42(17) 


6.42(20) 


B 5 


7.325 


10 


30 


5 


30.7(4) 


3.32(18) 


6.1(3) 


Ci 


7.234 


4 


16 


4 


27.6(3) 


57.44(25) 


7.18(3) 


c 2 


7.426 


5 


20 


5 


34.5(4) 


36.5(4) 


7.13(8) 


c 3 


7.584 


6 


24 


4 


41.4(5) 


24.7(4) 


6.94(12) 



Table 1: Lattice parameters and numerical results with z = (0,0,2a). 



If the "reweighting" observable 



is defined, then 



0[U,r i+1 }=e 



iT(z,r- 



S[U,U 4 ,r, +] ]-S[U,U 4 ,n] 



and the entropy can be computed as 



[0[U,r i+l ]) rj+i , 



2 V 1 , [ 3T(z,n) 

s = — — 2- ln 



z 2 r 2 L 3 



i=0 



3T(z,r i+1 ) 



(4.6) 



(4.7) 



(4.8) 



with z = (0,0,n z a) and with the integer n z being kept fixed when a — > 0. 

5. Results and conclusions 

Our goal is to show that the entropy density can be computed in the thermodynamic and con- 
tinuum limit by using Eq. (p~8|). To this aim we have calculated the entropy at three temperatures, 
1.5, 4.1 and 9.2r c , see Table |l| for the numerical results. The update algorithm used is a standard 
combination of heatbath and over-relaxation sweeps. The only changes over the standard algo- 
rithm reside (a) in the computation of the extra "staples" that determine the contribution to the 
action S[U,U4,r] of a given link variable U^{x), and (b) in the more frequent updating of the two 
time-slices on which the observable 0[t/,r,- + i] has its support. The bare coupling was tuned 



using the data of Ref. [1C] in order to match lattices of different Lo/a to the same temperature. 
Motivated by a study of the free case, we chose n z = 2. A comparison of the values of s for the 
lattices A\ and A\ a , and Bi and Bi fl indicates that finite- volume effects are very small or negligible 
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(s/T )=32tc745 



0.01 0.02 



0.03 0.04 
2„2 

a T 



• T=1.5T 
„ T=4.1 T 

♦ T=9.2T 



Figure 1: Scaling behavior of s/T , see Eq. (4.8). The Stefan-Boltzmann value is also displayed. 



within our statistical errors for these lattices. The behavior of the entropy density as a function of 
the lattice spacing is displayed in Figure [|. Cutoff effects are clearly quite mild. For illustration in 
the same plot we show a linear extrapolation in a 2 for the two lower temperatures where we have 
enough points. The corresponding intercepts are 4.77(8) and 6.30(9) at 1.5 and 4.1T C respectively, 
where the errors are statistical only. These results prove the numerical feasibility of the strategy 
discussed in the previous sections for computing the entropy density of the SU(3) Yang-Mills the- 
ory. It is also worth noting that, even though the statistical errors for the most expensive runs at 
Lq/q = 8, 10 are still quite large for a solid continuum extrapolation, the continuum-limit results at 



1.5 and 4. \T C are compatible with the best published ones Qll| , |12| ]. 

I thank Harvey B. Meyer for the intense and productive collaboration over the last year which 
led to the results discussed here. Many thanks to the organizers of the conference for their work, 
and for giving us the possibility to discuss physics in such a wonderful place. 
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